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Preface 


This  endeavor  was  initiated  to  determine  the  feasibility  of  fuel  analysis  using  Gas  Liquid-Phase 
Chromatography  (GLC).  The  GLC  is  one  of  several  analytical  techniques  that  is  under 
consideration  for  use  in  a  proposed  petroleum  product  analyzing  system. 
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SECTION  1 


INTRODUCTION 

An  earlier  report  detailed  the  use  of  GLC  as  an  effective  method  for  fuel  type  classification 
and  parameter  prediction.®1  In  that  study  sixty-four  fuel  samples  were  analyzed  using  a  single  GLC 
column  and  oven  parameter  profile  pair.  The  effectiveness  of  the  generated  models  at  fuel  type 
classification  and  fuel  parameter  prediction  were  established  by  statistical  analysis  of  the 
classification  and  prediction  errors. 

To  optimize  the  prediction  and  oven  profile  parameter  pairs,  seventeen  samples  were 
repetitively  analysis  using  seven  (7)  different  GLC  columns  and  oven  parameters.®  From  this  study 
the  DB-l/HP-1  columns  appeared  to  be  the  best  columns  for  the  current  analysis  as  was  evident  by 
the  minimum  prediction  errors  with  these  columns.  In  an  effort  to  reduce  the  analysis  time,  oven 
profiles  were  established  to  give  an  analysis  time  of  twelve  minutes.  The  prediction  errors  of  the 
model  generated  using  the  new  oven  profiles  were  compared  to  the  results  from  the  sixty-minute 
data  analysis. 


1  Superscript  numbers  ir  parenthesis  refer  to  references  at  the  end  of  the  report. 
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SECTION  2  TEST  PARAMETERS 


FUEL  SAMPLE  SET 

To  conduct  this  study,  sixty-seven  (67)  fuel  samples  were  used  of  which  fifteen  were  Jet  A-l 
or  Jet  A  conforming  to  ASTM  D1655,  "Standard  Specification  for  Aviation  Turbine  Fuels"(3), 
twelve  were  Grade  JP-4  and  nine  Grade  JP-5  fuels  conforming  to  Military  Specification  MIL-T- 
5624P,  "Turbine  Fuel,  Aviation,  Grades  JP-4,  JP-5,  and  JP-5/JP-8  ST"(4);  twelve  were  AVGAS 
conforming  to  ASTM  D  910,  "Standard  Specification  for  Aviation  Gasolines "(5);  nine  were  JP-8 
conforming  to  Military  Specification  MIL-T-83133,  "Turbine  Fuels,  Aviation,  Kerosene  Types, 
NATO  F-34  (JP-8)  and  NATO  F-35"<6);  seven  were  MOGAS  conforming  to  Military  Specification 
MIL-G-3056F,  "Gasoline,  Automotive,  Combat,,(7);  and,  three  Grade  DF-2  fuels  conforming  to 
Federal  Specification  VV-F-800D,  "Fuel  Oil,  Diesel"(8). 

The  fuel  samples  used  in  the  study  were  obtained  from  three  sources:  the  Fuels  Branch, 

Fuels  and  Lubricant  Division,  Wright-Patterson  Air  Force  Base  (WPAFB),  OH;  Directorate  of 
Aerospace  Fuels,  Detachment  13  (the  Kelly  AFB  contingent  at  WPAFB);  and  the  US  Army 
TARDEC  Fuels  and  Lubricants  Research  Facility  (TFLRF)  at  SwRI,  San  Antonio,  TX  (formerly  the 
Belvoir  Fuels  and  Lubricants  Research  Facility  (BFLRF)).  A  list  of  the  fuel  samples  used  in  the 
study  is  given  in  APPENDIX  A  and  the  specification  test  data  for  the  fuel  samples  is  given  in  the 
initial  report  (1). 


GLC  Column  and  Oven  Parameter  Pairs 

As  performed  in  previous  studies(1)(2),  the  fuel  samples’  GLC  chromatograms  were  reduced 
from  as  many  as  a  few  hundred  peaks  to  as  few  as  fifteen  peaks.  This  reduction  was  facilitated  by 
utilizing  a  peak-lumping  feature  in  the  chromatographic  collection  software.  The  peaks  were 
lumped  by  combining  all  of  the  peaks  in  a  given  time  interval  into  a  single  peak.  The  retention 
times  of  the  alkane  standards  and  the  time  intervals  for  the  chosen  regions  for  both  the  sixty-minute 
and  the  twelve-minute  analyses  are  given  in  Tables  1  and  2  respectively. 

CONDITIONS  OF  FID  Sixty-Minute  &  HP-1 

The  HP5890  GC  was  set  to  an  initial  temperature  of  40°C,  with  a  hold  time  for  5.0  minutes, 
a  program  ramp  6.0°C  per  minute,  to  a  final  temperature  of  250°C,  with  a  hold  time  of  20.0 
minutes.  A  Flame  Ionization  Detector  (FID)  was  used  with  a  30  meter  megabore  column,  0.53  cm 
inside  diameter  with  a  2.65  /*m  film  coat  of  methyl  silicone  (HP-1).  The  carrier  gas  was  helium 
with  a  flow  rate  of  8.6  mL  per  minute.  The  time  constraints  for  the  segmented  regions  and 
retention  times  for  the  alkane  standards  are  given  in  table  1: 
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TABLE  1:  Time  Constraints  For  FID  Sixty-Minute  &  HP-1  Analysis 


Region 

n-  Alkane 

Retention  Time 

Time  Boundaries 

1 

C5 

1.589  min 

0.00  -  2.27  min 

2 

C6 

2.937  min 

2.27  -  4.62  min 

3 

C7 

6.298  min 

4.62  -  8.33  min 

4 

C8 

10.368  min 

8.33  -  12.24  min 

5 

C9 

14.117  min 

12.24  -  15.77  min 

6 

CIO 

17.413  min 

15.77  -  18.95  min 

7 

Cll 

20.475  min 

18.95  -  21.87  min 

8 

C12 

23.261  min 

21.87  -  24.55  min 

9 

C13 

25.828  min 

24.55  -  27.07  min 

10 

C14 

28.315  min 

27.07  -  29.37  min 

11 

C15 

30.412  min 

29.37  -  31.54  min 

12 

C16 

32.663  min 

31.54  -  33.66  min 

13 

C17 

34.647  min 

33.66  -  35.64  min 

14 

C18 

36.645  min 

35.64  -  37.50  min 

15 

C19 

38.363  min 

37.50  -  39.26  min 

16 

C20 

40.152  min 

39.26  -  41.44  min 

17 

C21 

41.44  -  44.02  min 

18 

C22 

44.02  -  46.61  min 

19 

C23 

46.61  -  49.19  min 

20 

C24 

50.476  min 

49.19  -  60.00  min 

CONDITIONS  OF  FID  Twelve-Minute  &  HP-1 

The  HP5890  GC  was  set  to  initial  temperature  of  40°C,  with  a  hold  time  for  0.0  minutes,  a 
program  ramp  25.0°C  per  minute,  to  a  final  temperature  of  250°C,  with  a  hold  time  of  7.0  minutes 
A  FID  was  used  with  a  10  meter  megabore  column,  0.53  cm  inside  diameter  with  a  2.65  /xm  film 
coat  of  methyl  silicone  (HP-1).  The  carrier  gas  was  helium  with  a  flow  rate  of  9.6  mL  per  minute. 
The  time  constraints  for  the  segmented  regions  and  the  retention  times  for  the  alkane  standards  is 
given  in  table  2. 
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TAliLE  2:  Time  Constraints  For  FID  Twelve-Minute  &  HP-1  Analysis 


Region 

Retention  Time 

Time  Boundaries 

1 

C5 

0.591  min 

0.00  -  0.75  min 

2 

C6 

0.907  min 

0.75  -  1.16  min 

3 

Cl 

1.419  min 

1.16  -  1.74  min 

4 

C8 

2.068  min 

1.74 -  2.42  min 

5 

C9 

2.773  min 

2.42  -  3.12  min 

6 

CIO 

3.468  min 

3.12 -3.80  min 

7 

Cll 

4.136  min 

3.80  -  4.45  min 

8 

C12 

4.764  min 

4.45  -  5.06  min 

9 

C13 

5.355  min 

5.06  -  5.63  min 

10 

C14 

5.912  min 

5.63  -  6.17  min 

11 

C15 

6.434  min 

6.17  -  6.68  min 

12 

C16 

6.931  min 

6.68  -  7.16  min 

13 

C17 

7.397  min 

7.16  -  7.63  min 

14 

C18 

7.857  min 

7.63  -  8.07  min 

15 

C19 

8.284  min 

8.07  -  8.49  min 

! 

16 

C20 

8.693  min 

8.49  -  9.04  min 

17 

C21 

9.04  -  9.74  min 

18 

C22 

9.74  -  10.44  min 

19 

C23 

10.44  -  11.14  min 

20 

C24 

11.489  min 

11.14-  14.98  min 

SECTION  3  TEST  RESULTS  FOR  60-MINUTE  ANALYSIS 
Classification  Prediction 


From  the  initial  analysis,  samples  #  339  (MO-GAS),  93-F-591  (JP-4)  and  93-F-643 
(AVGAS)  were  identified  as  statistical  outliers  from  the  other  samples.  These  determinations  were 
based  on  a  Hierarchical  Cluster  Analysis  (HCA)  which  grouped  the  samples  based  on  the  similarity 
of  their  GLC  data.  With  a  HCA  analysis,  the  three  outliers  were  placed  into  grouping  outside  of 
the  other  samples  of  the  identical  type. 

Once  the  outliers  were  removed  HCA  analysis  resulted  in  five  groups  consisting  of  AVGAS, 
JP-4,  MO-GAS,  DF-2,  and  a  combined  group  of  JP-5,  JP-8,  and  JET  A.  There  were  no 
misclassifications  and  all  of  the  samples  were  assigned  to  their  proper  class.  The  class  assignments 
from  the  KNN  model  for  each  of  the  fuel  samples  are  given  in  Appendix  B. 

To  test  the  rigors  of  the  classification  model,  seventeen  randomly  selected  samples  were 
removed  from  the  data  set  and  a  KNN  model  was  made  using  the  remaining  data  points.  The 
seventeen  samples  included  eight  JP-8  group  samples,  three  JP-4,  three  AVGAS,  two  MOGAS,  and 
one  DF-2.  The  generated  classification  model  was  then  applied  to  the  removed  samples  and  resulted 
in  100  %  correct  class  determination  for  the  samples  in  the  removed  set.  A  list  of  the  removed 
samples  and  their  class  assignments  are  given  in  Appendix  C. 

PARAMETER  PREDICTION 

The  GLC  data  was  also  used  to  develop  models  to  predict  the  10  %,  50  %,  90  %,  and  final 
boiling  point  temperatures;  the  density;  and  the  viscosity  at  40°  and  100°C.  Using  the  results  of  the 
predictions,  the  average  error  and  average  %  error  were  calculated.  Table  3  is  a  listing  of  the 
average  error  and  average  %  error  for  the  generated  prediction  models. 


TABLE  3 :  Statistical  Analysis  of  Parameter  Predictions 


PARAMETER 

JP-8  SEPARATED  MODEL1 

TOTAL  MODEL2 

ERROR3 

%  ERROR4 

ERROR 

%  ERROR 

10%  BP  TEMP  (°K) 

5.14 

1.13 

10.2 

2.57 

50%  BP  TEMP  (°K) 

3.61 

0.74 

7.86 

1.80 

90%  BP  TEMP  (°K) 

6.39 

1.23 

11.7 

2.44 

FINAL  BP  TEMP  (°K) 

9.68 

1.77 

16.9 

3.29 

DENSITY  (g/mL) 

0.00440 

0.546 

0.00859 

1.10 

VIS  @  40°C  (cSt) 

0.191 

11.7 

0.235 

19.1 

VIS  @  100°C  (cSt) 

0.0973 

13.5 

0.0575 

10.4 

1 .  The  JP-8  separated  model  was  a  correlation  model  that  was  developed  from  samples  of  the  JP-8 
group. 

2.  The  total  model  was  a  correlation  model  that  was  developed  from  all  of  the  samples. 

N 

3.  The  average  error  =  where  Vcj  and  Vmi  are  the  calculated  and  the  measured  value 

of  the  Nth  parameter.  N 

( average  error)*  100 

4.  The  average  %  error  =  average  v~ 
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To  test  the  rigor  of  the  models  at  predicting  parameters  of  samples  that  weren’t  included  in 
the  correlation  models,  five  data  sets  were  removed  from  the  correlation  models  generation  sets  and 
Partial  Least  Squares  (PLS)  regression  method  were  used  to  predict  the  listed  parameters.  Tables  4 
thru  10  are  listing  of  the  predicted  values  and  errors  of  prediction  for  the  five  samples.  Plots  of  the 
predicted  versus  the  measured  parameters  for  all  samples  with  all  samples  included  in  the  correlation 
models  generation  sets  are  given  in  Appendix  D. 


TABLE  4:  10  %  Boiling  Point  Temperature 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(°K) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(°K) 

%  ERROR 

PREDICTED 

(°K) 

%  ERROR 

92-POSF-2959 

449 

439 

2.23 

443 

1.34 

92-POSF-2934 

424 

439 

3.53 

440 

3.77 

92-POSF-2747 

458 

459 

0.22 

468 

2.18 

93-F-304 

475 

456 

4.00 

449 

5.47 

93-F-668 

442 

458 

3.62 

456 

3.17 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 


TABLE  5.  50%  Boiling  Point  Temperature 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(°K) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(°K) 

%  ERROR 

PREDICTED 

(°K) 

%  ERROR 

92-POSF-2959 

472 

475 

0.63 

468 

0.85 

92-POSF-2934 

475 

478 

0.63 

465 

2.11 

92-POSF-2747 

466 

468 

0.43 

481 

3.22 

93-F-304 

505 

485 

3.96 

483 

4.36 

93-F-668 

484 

486 

0.41 

490 

1.24 

1.  The  actual  value  point  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  samples. 
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TABLE  6.  90%  Boiling  Point  Temperature 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(°K) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(°K) 

%  ERROR 

PREDICTED 

(°K) 

%  ERROR 

92-POSF-2959 

504 

516 

2.38 

511 

1.39 

92-POSF-2934 

522 

516 

1.15 

517 

0.96 

92-POSF-2747 

484 

481 

0.62 

483 

0.21 

93-F-304 

548 

533 

2.74 

524 

4.38 

93-F-668 

533 

536 

0.56 

523 

1.88 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 


Table  7.  Final  Boiling  Point  Temperature 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(°K) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(°K) 

%  ERROR 

PREDICTED 

(°K) 

%  ERROR 

92-POSF-2959 

523 

535 

2.29 

537 

2.68 

92-POSF-2934 

543 

527 

2.95 

543 

0.00 

92-POSF-2747 

509 

509 

0.00 

515 

1.18 

93-F-304 

566 

565 

0.18 

551 

2.65 

93-F-668 

600 

589 

1.83 

553  . 

7.83 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed 
using  data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  points  from  all  of  the  samples. 
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Table  8.  Density 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(g/mL) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(g/mL) 

%  ERROR 

PREDICTED 

(g/mL) 

%  ERROR 

92-POSF-2959 

792 

799 

0.90 

797 

0.63 

92-POSF-2934 

808 

796 

1.44 

796 

1.18 

92-POSF-2747 

807 

804 

0.36 

790 

2.11 

92-F-304 

812 

807 

0.62 

804 

0.99 

92-F-668 

823 

810 

1.58 

813 

1.22 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 


Table  9.  Viscosity  at  40° C 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(cSt) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(cSt) 

%  ERROR 

PREDICTED 

(cSt) 

%  ERROR 

92-POSF-2959 

1.442 

1.510  1 

4.716 

1.432 

0.694 

92-POSF-2934 

1.545 

1.498 

3.042 

1.591 

2.977 

92-POSF-2747 

1.530 

1.526 

0.261 

1.244 

18.69 

93-F-304 

1.769 

1.715 

3.053 

1.926 

8.875 

93-F-668 

1.897 

1.725 

9.067 

1.878 

1.002 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 
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Table  10.  Viscosity  at  100°C 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(cST) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(cSt) 

%  ERROR 

PREDICTED 

(cSt) 

%  ERROR 

92-POSF-2959 

0.7922 

0.7914 

0.101 

0.7411 

6.450 

92-POSF-2934 

0.7586 

0.7420 

2.188 

0.7296 

3.823 

92-POSF-2747 

0.7464 

0.6540 

12.38 

0.7795 

4.435 

93-F-304 

0.8248 

0.9122 

10.60 

0.8289 

0.497 

93-F-668 

0.8169 

0.9346 

14.41 

0.9051 

11.08 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 
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SECTION  4  TEST  RESULTS  FOR  TWELVE-MINUTE  ANALYSIS 
CLASSIFICATION  PREDICTION 

From  the  initial  analysis,  sample  #  93-F-643  (AVGAS)  was  identified  as  statistical  outlier 
from  the  other  samples.  This  determination  was  based  on  a  Hierarchical  Cluster  Analysis  (HCA) 
which  grouped  the  samples  based  on  the  similarity  of  their  GLC  data.  With  a  HCA  analysis, 
sample  #  93-F-643  was  placed  into  a  group  outside  the  other  samples  of  identical  type. 

Once  the  outlier  was  removed,  HCA  analysis  resulted  in  four  groups  consisting  of  AVGAS, 
DF-2,  a  combined  group  of  JP-4  and  MOGAS,  and  a  combined  group  of  JP-5,  JP-8,  and  JET  A. 
Three  samples  were  misclassifications  with  all  of  the  other  samples  being  to  their  proper  class.  It 
should  be  noted  that  for  each  samples  misclassified,  a  duplicate  GLC  analysis  was  conducted  and 
was  properly  classified.  The  class  assignments  from  the  KNN  model  for  each  of  the  fuel  samples 
are  given  in  Appendix  B. 

To  test  the  rigors  of  the  classification  model,  seventeen  randomly  selected  samples  were 
removed  from  the  data  set  and  a  KNN  model  was  made  using  the  remaining  data  points.  The 
seventeen  samples  included  eight  JP-8  group  samples,  five  JP-4  group  samples,  three  AVGAS,  and 
one  DF-2.  The  generated  classification  model  was  then  applied  to  the  removed  samples  and  resulted 
in  all  of  the  fuel  samples  being  correctly  classified.  A  list  of  the  removed  samples  and  their  class 
assignments  are  given  in  Appendix  C. 


PARAMETER  PREDICTION 

The  GLC  data  was  also  used  to  develop  models  to  predict  the  10  %,  50  %,  90  %,  and  final 
boiling  point  temperatures;  the  density;  and  the  viscosity  at  40°  and  100°C.  Using  the  results  of  the 
predictions,  the  errors  were  calculated.  Table  1 1  is  a  listing  of  the  average  error  and  average  % 
error  for  the  generated  models. 


TABLE  1 1 :  Statistical  Analysis  of  Parameter  Predictions 


PARAMETER 

.TP-8  SEPARATED  MODEL1 

TOTAL  M< 

ODEL2 

ERROR3 4 

ERROR 

%  ERROR 

10%  BP  TEMP  UK) 

6.19 

1.37 

12.0 

3.23 

50%  BP  TEMP  f°K) 

5.33 

1.10 

9.82 

2.44 

90%  BP  TEMP  (°K) 

9.47 

1.85 

14.9 

3.23 

FINAL  BP  TEMP  (°K) 

13.3 

2.42 

20.3 

4.01 

DENSITY  (s/mL) 

0.00361 

0.45 

0.00596 

0.76 

0.198 

11.8 

0.218 

17.1 

0-0512 _ 

_ L55 _ _ _ 

IHRSSSi 

_ Zte, _ _ 

1 .  The  JP-8  separated  moc 

lei  was  a  correlation  model  that  was  developed  1 

rom  samples  of 

group. 

2.  The  total  model  was  a  correlation  model  that  was  developed  from  all  of  the  samples. 


3.  The  average  error  =  P  Vd~Vmi\  where  Vci  and  Vmi  are  the  calculated  and  the  measured  value  of 

the  Nth  parameter.  V 

(i average  error)  *100 

4.  The  average  %  error  =  average  Vmi 
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To  test  the  rigor  of  the  models  at  predicting  parameters  of  samples  that  weren’t  included  in 
the  correlation  model,  five  data  sets  were  removed  from  the  correlation  models  and  Partial  Least 
Squares  (PLS)  regression  method  were  used  to  predict  the  listed  parameters.  Tables  12  thru  18  are 
listing  of  the  predicted  values  and  errors  of  prediction  for  the  five  samples.  Plots  of  the  predicted 
versus  the  measured  parameters  for  all  samples  with  all  samples  included  in  the  correlation  model 
are  given  in  Appendix  D. 


TABLE  12:  10  %  Boiling  Point  Temperature 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(°K) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(°K) 

%  ERROR 

PREDICTED 

(°K) 

%  ERROR 

92-POSF-2959 

449 

450 

0.22 

435 

3.12 

92-POSF-2934 

424 

451 

6.37 

434 

2.36 

92-POSF-2747 

458 

444 

3.06 

425 

7.21 

93-F-304 

475 

456 

4.00 

449 

5.47 

93-F-668 

442 

451 

2.04 

456 

3.17 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 


TABLE  13.  50%  Boiling  Point  Temperature 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(°K) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(°K) 

%  ERROR 

PREDICTED 

(°K) 

%  ERROR 

92-POSF-2959 

472 

482 

2.12 

475 

0.64 

92-POSF-2934 

475 

481 

1.26 

471 

0.84 

92-POSF-2747 

466 

488 

4.72 

511 

9.66 

93-F-304 

505 

484 

4.16 

487 

3.56 

93-F-668 

484, 

485 

0.21 

485 

0.21 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 
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TABLE  14.  90%  Boiling  Point  Temperature 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(°K) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(°K) 

%  ERROR 

PREDICTED 

(°K) 

%  ERROR 

92-POSF-2959 

504 

521 

3.37 

532 

5.56 

92-POSF-2934 

522 

521 

0.19 

531 

1.72 

92-POSF-2747 

484 

524 

8.26 

551 

13.8 

93-F-304 

548 

520 

5.11 

524 

4.38 

93-F-668 

533 

521 

2.25 

530 

0.56 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  Tne  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 


Table  15.  Final  Boiling  Point  Temperature 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(°K) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(°K) 

%  ERROR 

PREDICTED 

(°K) 

%  ERROR 

92-POSF-2959 

523 

545 

2.29 

567 

8.41 

92-POSF-2934 

543 

542 

0.20 

566 

4.24 

92-POSF-2747 

509 

562 

10.4 

593 

16.5 

93-F-304 

566 

536 

5.30 

547 

3.36 

93-F-668 

600 

552 

8.00 

563 

6.17 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 
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Table  16.  Density 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(g/mL) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(g/mL) 

%  ERROR 

PREDICTED 

(g/mL) 

%  ERROR 

92-POSF-2959 

792 

806 

1.77 

805 

1.64 

92-POSF-2934 

808 

806 

0.25 

805 

0.37 

92-POSF-2747 

807 

809 

0.37 

813 

0.74 

92-F-304 

812 

806 

0.74 

807 

0.62 

92-F-668 

823 

807 

1.94 

808 

1.82 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 


Table  17.  Viscosity  at  40° C 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(cSt) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(cSt) 

%  ERROR 

PREDICTED 

(cSt) 

%  ERROR 

92-POSF-2959 

1.442 

1.703 

18.10 

1.815 

25.87 

92-POSF-2934 

1.545 

1.716 

11.07 

1.792 

15.99 

92-POSF-2747 

1.530 

1.535 

0.327 

1.518 

0.784 

93-F-304 

1.769 

1.780 

0.622 

1.679 

5.088 

93-F-668 

1.897 

1.702 

10.28 

1.674 

11.76 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 
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Table  18.  Viscosity  at  100°C 


SAMPLE 

NUMBER 

ACTUAL 

VALUE1 

(cSt) 

GROUP  1  MODEL2 

TOTAL  MODEL3 

PREDICTED 

(cSt) 

%  ERROR 

PREDICTED 

(cSt) 

%  ERROR 

92-POSF-2959 

0.7922 

0.7584 

4.267 

0.8024 

1.288 

92-POSF-2934 

0.7586 

0.7616 

0.395 

0.7932 

4.561 

92-POSF-2747 

0.7464 

0.7187 

3.711 

0.7237 

3.041 

93-F-304 

0.8248 

0.8075 

2.0975 

0.8156 

1.115 

93-F-668 

0.8169 

0.7802 

4.493 

0.7775 

4.823 

1.  The  actual  value  is  the  parameter’s  measured  value. 

2.  The  Group  1  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using 
data  point  from  the  JP-8  group  samples. 

3.  The  Total  Model  is  the  parameter’s  predicted  value  using  a  model  that  was  developed  using  data 
points  from  all  of  the  samples. 
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SECTION  V 
CONCLUSION 


For  both  the  twelve-minute  and  the  sixty-minute  analysis,  two  set  of  prediction  were 
determined.  One  set  of  prediction  data  consisted  of  using  models  that  were  generated  from  JET  A, 
JET  A-l ,  JP-5  and  JP-8  (or  collectively  termed  Group  1  samples)  samples  and  the  second  set  of 
prediction  data  consisted  of  using  models  that  were  generated  from  all  sixty-seven  fuel  samples. 

With  both  the  twelve-minute  and  the  sixty-minute  prediction  set,  the  average  percent  errors 
of  the  predicted  versus  the  actual  parameter  values  are  relatively  small  (less  than  2  %).  Also,  the 
difference  between  the  average  percent  errors  of  model  developed  from  GROUP  1  samples  and  the 
models  developed  from  all  sixty-seven  data  sets  were  relative  small  (less  than  2  %)  for  the  10  %,  50 
%,  90  %,  and  final  boiling  point  temperature;  and  the  density.  The  low  errors  for  both  the  sixty- 
minute  and  twelve-minute  analyses  suggest  that  both  model  sets  are  useful  for  the  prediction  of  the 
above  parameters. 

Although  the  errors  for  the  twelve-minute  analyses  are  slightly  higher  than  the  errors  of  the 
sixty-minute  analyses,  the  reduction  in  sample  analysis  time  more  than  makes  up  for  the  slightly 
higher  errors.  This  shows  that  twelve-minute  analysis  time  may  be  used  to  decrease  sample  test 
time  and  to  increase  sample  through  put. 

For  the  predictions  of  the  viscosities  at  40°  and  100°C  for  both  analyses,  the  average  percent 
error  of  the  predicted  versus  the  actual  parameter  values  are  relatively  large  (more  than  10  %) 
which  shows  that  these  parameters  are  not  good  candidates  for  prediction  using  the  GLC  data. 

Future  work  will  involve  studing  the  effects  of  column  wall  thickness  on  the  GLC  analysis, 
and  inclusion  of  additional  fuel  samples  and  parameters  for  prediction. 
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APPENDIX  A 


List  of  Fuel  Samples 


A1 


A2 


TABLE  19.  List  of  Fuel  Samples 


SAMPLE  ORIGIN 


WPAFB  BLDG  490 
WPAFB  BLDG  490 
WPAFB  BLDG  490 
WPAFB  BLDG  490 
WPAFB  BLDG  490 
WPAFB  BLDG  490 
WPAFB  BLDG  490 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  BLDG  490 
WPAFB  BLDG  490 
WPAFB  BLDG  490 
WPAFB  BLDG  490 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 
WPAFB  KELLY  AFB  DET 


FUEL  TYPE 

TEST  DATA 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A 

YES 

JET  A-l 

YES 

JP-8 

YES 

JP-8 

YES 

JP-5 

YES 

JP-5 

YES 

JP-5 

YES 

JP-5 

YES 

JP-5 

YES 

JP-5 

YES 

JP-5 

YES 

JP-5 

YES 

DF-2 

YES 

JP-4 

YES 

JP-4 

YES 

JP-4 

YES 

JP-4 

YES 

JP-4 

YES 

JP-4 

YES 

JP-4 

YES 

JP-4 

YES 

JP-4 

YES 

JP-4 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

AVGAS 

YES 

MOGAS 

YES 

MOGAS 

YES 

MOGAS 

YES 

MOGAS 

YES 

MOGAS 

YES 

MOGAS 

YES 

MOGAS 

YES 
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Classification  Prediction 


B1 


B2 


Table  20:  Result  of  Classification  Analysis  Using  12  &  60  min  data. 


SAMPLE  # 


93-POSF-2959 


92-POSF-2928 


92-POSF-2926 


93-POSF-2747 


92-POSF-2930 


91-POSF-2827 


92-POSF-2922 


93-F-173 


93-F-142 


93-F-304 


93-F-280 


93-F-665 


93-F-444 


93-F-560 


92-POSF-2931 


AL-19850-F 


AL-20123-F 


AL-20336-F 


AL-19903-F 


93-F-351 


92-POSF-2934 


AL-20335-F 


AL-20011-F 


92-POSF-2936 


93-POSF-2963 


91-POSF-2817 


93-F-284 


93-F-311 


FUEL 

TYPE 

Classification  of  samples  using  the  60 
min  data  analysis 

Classification  of  samples  using  the 

12  min  data  analysis 

JET  A 

GROUP  l1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A 

GROUP  1 

GROUP  1 

JET  A-l 

_ 

GROUP  1 

GROUP  1 

GROUP  1 

GROUP  1 

GROUP  1 

GROUP  1 

JP-8 

GROUP  1 

GROUP  1 

JP-8 

GROUP  1 

GROUP  1 

mm 

GROUP  1 

GROUP  1 

GROUP  1 

GROUP  1 

JP-8 

GROUP  1 

GROUP  1 

JP-8 

GROUP  1 

GROUP  1 

JP-8 

GROUP  1 

GROUP  1 

JP-5 

GROUP  1 

GROUP  1 

JP-5 

GROUP  1 

GROUP  1 

JP-5 

GROUP  1 

GROUP  1 

JP-5 

GROUP  1 

GROUP  1 

1  GROUP  1  consists  of  JP-5,  JP-8,  JET  A,  and  JET  A-l 
samples. 
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TABLE  20:  Result  of  Classification  Analysis  Using  12  &  60  min  data  (continued). 


SAMPLE#  I  FUEL  TYPE 


93-F-312 


93-F-310 


93-F-313 


93-F-374 


93-F-668 


93-F-201 


93-F-295 


93-F-289 


93-F-290 


93-F-279 


93-F-643 


93-F-551 


93-F-644 


93-F-338 


93-F-609 


93-F-610 


93-F-539 


93-L-100 


AL-20221-F 


AL-19915-F 


93-F-152 


DEJESS 


93-F-347 


93-F-412 


93-F-402 


93-F-233 


93-F-586 


93-F-624 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


DF-2 


DF-2 


DF-2 


Classification  of  samples  using 
the  60  min  data  analysis 


GROUP  1 


GROUP  1 


GROUP  1 


GROUP  1 


GROUP  1 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


DF-2 


DF-2 


Classification  of  samples  using 
the  12  min  data  analysis 


GROUP  1 


GROUP  1 


GROUP  1 


GROUP  1 


GROUP  1 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


AVGAS 


DF-2 


DF-2 


DF-2 


GROUP  42 


GROUP  4 


GROUP  4 


GROUP  4 


GROUP  4 


GROUP  4 


GROUP  4 


GROUP  4 


2  GROUP  4  consists  of  JP-4  and  MOGAS 


B4 


TABLE  20:  RESULT  OF  CLASSIFICATION  ANALYSIS  USING  12  &  60  MIN 
DATA  (continued). 


SAMPLE# 

FUEL 

TYPE 

Classification  of  samples  using 
the  60  min  data  analysis 

Classification  of  samples 
using  the  12  min  data 
analysis 

93-F-625 

JP-4 

JP-4 

GROUP  4 

93-F-640 

JP-4 

JP-4 

GROUP  4 

93-F-591 

JP-4 

JP-4 

GROUP  4 

93-F-653 

JP-4 

JP-4 

GROUP  4 

93-F-326 

MOGAS 

MOGAS 

GROUP  4 

93-F-307 

MOGAS 

MOGAS 

GROUP  4 

93-F-339 

MOGAS 

MOGAS 

GROUP  4 

93-F-306 

MOGAS 

MOGAS 

GROUP  4 

93-F-464 

MOGAS 

MOGAS 

GROUP  4 

93-F-638 

MOGAS 

MOGAS 

GROUP  4 

93-F-449 

MOGAS 

MOGAS 

GROUP  4 
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Classification  Prediction  of  Removed  Samples 


c-i 


C-2 


Table  21:  Classification  of  Removed  Samples  Using  12  &  60  min  data 


SAMPLE# 

FUEL  TYPE 

Classification  of  samples 
using  the  60  min  data 
analysis 

Classification  of  samples 
using  the  12  min  data 
analysis 

92-POSF-2928 

JET  A 

GROUP  l3 

GROUP  1 

93-F-444 

JET  A 

GROUP  1 

GROUP  1 

93-F-665 

JET  A 

GROUP  1 

GROUP  1 

93-F-280 

JET  A 

GROUP  1 

GROUP  1 

92-POSF-2934 

JP-8 

GROUP  1 

GROUP  1 

AL-20335-F 

JP-8 

GROUP  1 

GROUP  1 

93-F-311 

JP-5 

GROUP  1 

GROUP  1 

93-F-310 

JP-5 

GROUP  1 

GROUP  1 

93-F-289 

AVGAS 

AVGAS 

AVGAS 

93-F-290 

AVGAS 

AVGAS 

AVGAS 

93-F-201 

AVGAS 

AVGAS 

AVGAS 

93  F-640 

JP-4 

JP-4 

GROUP  44 

DEJESS 

JP-4 

JP-4 

GROUP  4 

93-F-402 

JP-4 

JP-4 

GROUP  4 

93-F-307 

MOGAS 

MOGAS 

GROUP  4 

93-F-326 

MOGAS 

MOGAS 

GROUP  4 

AL-20221-F 

DF-2 

DF-2 

DF-2 

3  GROUP  1  consists  of  JP-5f  JP-8,  JET  A,  and  JET  A-l 
samples. 

4  GROUP  4  consists  of  JP-4  and  MOGAS  samples. 
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APPENDIX  D 


Plots  of  Measured  Versus  Predicted  Parameters 


D1 


10%  BP  TEMP  (60  MINUTE  ANALYSIS) 

R2=0.9494 


D2 


PREDICTED  TEMP.  (°K) 


600 


90%  BP  TEMP  (60  MINUTE  ANALYSIS) 

R=0.9326 


MEASURED  TEMP.  (°K) 


D 


PREDICTED  DENSITY  (g/mL) 


DENSITY  (60  MINUTE  ANALYSIS) 

r!=0.9446 


MEASURED  DENSITY  (g/mL) 


D6 


VISCOSITY  @  100°C  (60  MINUTE  ANALYSIS) 

R2=0.7163 


MEASURED  VISCOSITY  (cSt) 


D8 


PREDICTED  TEMP.  (°K) 


50%  BP  TEMP.  (12  MINUTE  ANALYSIS) 

R2=0.9143 


MEASURED  TEMP.  (°K) 


DIO 


PREDICTED  TEMP.  (°K) 


FINAL  BP  TEMP.  (12  MINUTE  ANALYSIS) 

R2=0.8014 


MEASURED  TEMP.  (°K) 


D12 
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